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Abstract. The morphological properties of the Carina, Sculptor and Fornax dwarfs are investigated using new 
wide field data with a total area of 29 square degrees. The stellar density maps are derived, hinting that Sculptor 
possesses tidal tails indicating interaction with the Milky Way. Contrary to previous studies we cannot find any 
sign of breaks in the density profiles for the Carina and Fornax dwarfs. The possible existence of tidal tails in 
Sculptor and of King limiting radii in Fornax and Carina are used to derive global M/L ratios, without using 
kinematic data. By matching those M/L ratios to kinematically derived values we are able to constrain the orbital 
parameters of the three dwarfs. Fornax cannot have M/L smaller than 3 and must be close to its perigalacticon 
now. The other extreme is Sculptor that needs to be on an orbit with an eccentricity bigger than 0.5 to be able 
to form tidal tails despite its kinematic M/L. 



1. Introduction 

Dwarf spheroidal galaxies (dSphs) are generally charac- 
terized by low masses (smaller than 10 s Mq), low surface 
brightnesses, low central concentrations, old populations, 
and the absence of interstellar gas. Due to their intrinsic 
faintness, their study has been and still is very difficult. 

Despite these apparent "signs of insignificance" , dSphs 
play an important role in Cold Dark Matter (CDM) sce- 
narios, where larger galaxies are built out of smaller build- 
ing blocks. Comparatively high velocity dispersions are 
found in dSphs which could be the signature of a fairly 
large amount of dark matter inside the visible radius. 
Indeed, mass to light (M/L) ratios derived from virial 
equilibrium can be as high as 100 for some dSphs (Mateo 
1998). Dwarf spheroidals would thus represent the low 
mass end of the mass function of dark matter halos. Yet 
the observed number of dwarf galaxies is by two orders of 
magnitudes too small compared to the prediction of Cold 
Dark Matter simulations (Klypin et al. 119991 Moore et al. 
1999}. This is a major challenge for CDM scenarios, which 
are very successful in explaining the large scale structure 
of the universe. There have been several propositions to 
alleviate this problem, which mainly deal with feedback 
effects or changes to CDM theory. Recently, D'Onghia et 
al. 1|2002[1 have shown that self-interacting dark matter 
fails to solve this abundance discrepancy. More conser- 
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vative solutions include feedback and heating by an ion- 
izing UV background, thereby separating the dark from 
the luminous matter. Along this line of thought Stoehr et 
al. I|2002[) identified the Galactic satellites with the most 
massive satellite substructures in their CDM simulations 
and inferred that tidal tails and extra-tidal stars should 
not, then be present in most systems because they are 
embedded in larger dark matter halos. 

On the other hand, dwarf galaxies without dark mat- 
ter can be produced naturally in the course of mergers be- 
tween bigger galaxies (Barnes Ll99„2j); they form as clumps 
in tidal tails. In the absence of dark matter, the Galactic 
dSphs will be disturbed by the tides from the gravita- 
tional potential of the Milky Way. Their large stellar ve- 
locity dispersions would then either be the result of the 
line-of-sight extension of the galaxy or denote unbound 
stars following the same orbits for a number of orbital pe- 
riods after disruption of the dwarf galaxy. Then the mor- 
phology of the dSphs should be disturbed and tidal tails 
should be visible. Kroupa l|1997f) and Klessen & Kroupa 
( 199SJ) have carried out simulations of dwarf galaxies with- 
out dark matter in the potential of a larger parent galaxy. 
They are able to produce galaxies with similar features as 
those observed for some dSphs. A conclusive test to distin- 
guish between the tidal and the dark matter model could 
come from analysis of the horizontal branch morphology 
of stars of the dSph galaxy in the color magnitude diagram 
(Klessen & Zhao l27)0"2*f . To what extend the presence of 
tidal tails would rule out dark matter is not entirely clear. 
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Mayer et al. ( 2001 ) state that, given a flat core for the DM 
halo, the production of tidal tails is possible. Yet Burkert 
(1996) argued that "The existence of extratidal stars in 
Carina, Draco and Ursa Minor would demonstrate that 
these dSphs cannot contain significant amounts of dark 
matter." 

Therefore it is of great interest to investigate the 
morphology of the satellites of the Milky Way. Irwin 
and Hatzidimitriou (1995, hereafter IH) used star count- 
ing techniques to investigate the morphology of all eight 
Galactic dSphs known at that time. They found that their 
stellar density profiles are well fitted by an empirical King 
profile (King 1962) in the inner parts. Yet they find for al- 
most all dSphs a departure from the same profile (or break 
in the profile) at several core radii, indicating a "break 
population" of stars. As the King profile is thought to 
be the density profile for a relaxed satellite in an external 
potential, this has often been taken as evidence for the ex- 
istence of tidal tails or extended components of unknown 
origin. 

Recently, Odenkirchen et al. (2001) have used data 
from the Sloan Digital Sky Survey to investigate the mor- 
phology of the Draco dwarf. Due to the immense sky cov- 
erage and the availability of five photometric bands, they 
were able to carry out a more sophisticated analysis. They 
defined a region in colour-magnitude space holding chiefly 
Draco stars. By selecting only the stars in this region, 
they were able to reduce the "background" of objects - 
consisting of foreground stars and background galaxies - 
by one order of magnitude. They could not confirm the 
existence of extratidal stars claimed in IH and Piatek et 
al. H2001(l . they see Draco as a fully relaxed object down to 
their detection limit. Tidal models for Draco can conclu- 
sively be ruled out on basis of the small dispersion of its 
horizontal-branch stars in the color-magnitude diagram 
(Klessen, Grebel, & Harbeck 2003). A similar technique 
was applied to the Ursa Minor dSph by Martinez-Delgado 
et al. H2UU1|I . They found clear indications for tidal inter- 
action between this dwarf and the Milky Way. Majewski 
et al. used even stronger cuts in colour-magnitude 

space to identify Red Giant Branch stars belonging di- 
rectly to each galaxy. They studied seven objects, five 
dSphs and two globular clusters and found a break in the 
radial RGB number density profile for every single object. 

This paper presents new wide and deep imaging of 
three dSphs, namely Carina, Sculptor and Fornax, aimed 
at searching for extratidal stars. The data is presented in 
section |2 In section we describe the derived isoplcth 
plots and morphological parameters for all three galaxies. 
In section 0] we discuss the results and draw our conclu- 
sions. 

2. Observations and reduction 

We obtained V band images of the three dwarf spheroidal 
galaxies in the constellations of Carina, Sculptor and 
Fornax with the Wide Field Imager (WFI) on the 
MPG/ESO 2.2m telescope on La Silla. The observing run 



took place on eleven nights between the 17th and the 29th 
of September 1999. The WFI provides a mosaic of 8 CCDs 
with a total FOV of 32' x 32'. The filling factor is 96% due 
to small gaps between the chips. We covered 4 square de- 
grees for Carina, 8.5 square degrees for Fornax and 16.25 
square degrees for Sculptor. To cover these large fields, 
we made no attempt to dither, so that the gaps between 
the WFI chips show up in the subsequent analysis. The 
typical exposure time per field is 600s, reaching a limiting 
magnitude of 23.5 mag on photometric nights. We also do 
have some I and B band images which are not complete 
enough to be used, but will be supplemented hopefully 
next September. 

The frames were bias-subtracted and flatfielded in the 
usual way and we applied a median filter for cosmic ray re- 
moval. Between 2 and 6 frames were taken per field, either 
one after another or on different nights. For each field they 
were aligned and coadded to reach fainter limiting mag- 
nitudes. The object catalogue was extracted from these 
coadded frames using the SExtractor software of E. Bertin 
(1996). The fluxes were determined in elliptical apertures 
around the center of each object. As the precision of the 
photometry is not critical for this work, we considered it 
unnecessary to spend the large amount of time required 
to do PSF-photometry. 

The atmospheric extinction was measured from 
Landolt (1992) standard star fields. Extinction coefficients 
where determined on a night to night basis, but as some 
object fields were observed two times on different nights we 
applied only a medium extinction correction to all frames. 
The mean extinction coefficient is 0.137 ±0.019. The final 
error in absolute photometry for the faintest stars is about 
0.3 mag, while the brighter stars (below 20th magnitude) 
are precise to about 0.1 mag. No correction fo possible 
varying galactic extinction over the fields has been made. 

Using the USNO-A1.0 catalogue we determined the 
astrometric calibration for each coadded frame. Then the 
object catalogues from each field were combined to a fi- 
nal catalogue. We used the overlapping parts of the fields 
(~ 15 arcmin 2 ) to correct for photometric offsets be- 
tween the fields. Due to the long duration of the observ- 
ing run, conditions varied from photometric to thin cirrus. 
Nonetheless typical offset values are lower than 0.1 mag. 
The way we deal with non-uniform limiting magnitudes 
for the isopleth maps and the density profile determina- 
tion is described in the relevant paragraphs. 

SExtractor also produces a classification parameter 
for the detected objects which was used to discriminate 
between galaxies and stars and thus reduce the number 
of background objects. This parameter runs between 
(galaxies) and 1 (stars). While there is a very clear di- 
chotomy above 21st magnitude, stars and galaxies are al- 
most indistinguishable below the limiting magnitude in 
each frame. From the distribution of the classification pa- 
rameter and eyeball inspection we adopted a cut at 0.32 
uniformly for all frames. This excludes between 4% and 
10% of the objects, depending on how many stars the dSph 
itself contributes to the whole count. We verified that the 



Walcher et al.: Morphology and dark matter in dSphs 



3 



excluded objects are distributed uniformly over the Carina 
field. Unfortunately this is not true for the Sculptor and 
Fornax fields as the seeing varied between 1" and 2.5" over 
the run and the discrimination is naturally hampered by 
bad seeing. We therefore choose to include all objects in 
the Sculptor and Fornax isopleth maps. For a discussion 
of the relevant effects see next chapter. 

3. Derivation of morphological parameters and 
determination of the profile 

From the final object catalogues we produced the stellar 
density maps shown as isopleth plots in Figures 1112*1 andOl 
To reduce the effects of the gaps between the WFI chips 
we first produced higher resolution density maps, where 
the gaps can be easily identified as vertical or horizontal 
rows with a density value of zero. We filled in these parts 
by taking the simple average of the surrounding pixels, ex- 
cluding the ones already defined as part of the gaps. We 
then rebinned all density maps to have a pixel size of f 
arcmin . This procedure does not add any information to 
the images but prevents the eye being mislead by spurious 
artefacts. While we used 23.5 mag as the flux limit, some 
frames have a shallower limiting magnitude than others 
due to differing observing conditions. Their locations on 
the maps are pointed out in the text. As their shape is 
rectangular, they can hardly be confused with real struc- 
ture (see for example the lower left corner of the Carina 
map). 

The final object catalogues contain foreground 
Galactic stars, stars from the dwarf galaxy itself and back- 
ground galaxies. These three components all have different 
distributions in apparent magnitude. We took advantage 
of this additional information by subdividing the object 
catalogue into magnitude bins, indexed by i in the follow- 
ing. We then produced density maps imgi in each mag- 
nitude bin. We also produced apparent magnitude his- 
tograms of a major part of the background area (called 
backi) and one area in the galaxy center (gali). We then 
determined weights in the following way: 



In coadding the density maps imgi to the final density 
map imgfin we used the following formula: 

imgfin = (imgi — backi) x wf al (2) 

This improves the contrast in surface density between 
the center of the galaxy and the background area by ~ 
20%. 

The density values of the contour lines in Figs, [l] to 
U3 have been chosen to be the background value (dotted 
line), 1(7 in the sense of statistics on a single pixel above 
that (thin solid line), 2a, 5a, 10a and so on (thick solid 
lines). The images have been smoothed to clear up the 
plots. 



For Carina we used a particuarly low smoothing fac- 
tor of 2 to be sure not to hide any signs of irregularity. 
We applied the star/galaxy classification. No significant 
departure from the spheroidal shape can be seen. In the 
southeast corner there is one shallower frame. A galactic 
gradient can be seen from the northeastern to the south- 
western corner. 

The Sculptor map is more difficult to interpret, be- 
cause of the varying seeing conditions over the run. For 
consistency over the field we show the density map pro- 
duced without star/galaxy discrimination. Sculptor looks 
very regular in the inner parts, with almost zero elliptic- 
ity. IH already noted an ellipticity increase with radius. On 
the outskirts two opposite annexes can be seen, where the 
eastern one obviously extends further to the south. These 
are hints but no proofs of tidal extensions, because the 
eastern annex actually disappears when we produce the 
same map taking into account the classification cut. This 
could either be, because the exceptional seeing on this field 
allows the classifier to identify more galaxies or because 
of a background cluster of galaxies. The reverse is true for 
the other annex where the seeing was exceptionally bad. 
This potential tail lies on a field with particularly bad see- 
ing. Also the fields surrounding this interesting part of the 
density map have somewhat shallower limiting magnitude 
due to weather, so that we can not follow the potential 
tails further out. We proposed further observations to ob- 
tain a second color which will lift this uncertainty. 

The same problem of seeing makes it uncertain to ap- 
ply the classification also for Fornax. The differences be- 
tween the maps are negligible, we therefore choose to show 
the map that includes the galaxies to avoid inconsistency. 
Fornax' innermost isopleths show a steeper decline to the 
east than to the west, as noticed earlier by Demers et al. 
(1994). Fornax is almost as extended as the whole field 
as shown by Fig 7. As can also be seen from the straight 
line in the isophotes, the frame south of the central field is 
shallower, as well as two further frames in the southeast- 
ern corner. These frames were masked out in the density 
profile derivation below. The apparent overdensity at the 
lower 39.6 right ascencsion tickmark is an edge effect of 
the gap interpolation routine. 

Centers, position angles and ellipticities were derived 
on the isopleth maps using the algorithm published in 
Bender & Mollenhof and are listed in Tabled They 

are generally in good agreement with the values published 
in IH. It is worth noting that the ellipticity of Sculptor in- 
creases with radius, from 0.15 in the very center to 0.3 at 
a radius of about 30 arcmin. 

To determine the density profile we went back to the 
object catalogues. We counted all objects brighter than 
23.5mag on logarithmically spaced concentric annuli of 
fixed position angle and ellipticity. We did not vary these 
parameters over radius. The background density was de- 
termined from those parts of the profile where the density 
reaches a constant value and then subtracted from the 
profile. We checked this value on a background area in 
the density map, both agree. We did not use the weight- 
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Fig. 1. Contour plot of the Carina dwarf spheroidal. The 
density levels correspond to background value (dotted 
line), la above that (thin solid line), 2a, 5a, 10<r and so 
on (thick solid lines). No significant departure from the 
spheroidal shape can be seen. A galactic gradient can be 
seen from the northeastern to the southwestern corner. 



ing scheme used in the production of the density maps 
but did use the star/galaxy classifier for the Carina pro- 
file. We masked out the frames with a shallower limiting 
magnitude discussed earlier. 

Though crowding is not dominant in our CCD images, 
we determined and applied a correction adapted to the ob- 
ject extraction mechanism of SExtractor. By construction 
the smallest separation where the detection algorithm of 
SExtractor resolves two stars is 1 FWHM. Therefore one 
star covers an area of it ■ FWHM 2 in the frame where no 
other star can be detected. The area covered by stars in 
arcmin 2 is therefore 



A, 



it ■ FWHM ■ n*, 



(3) 



where n* is the number of stars per arcmin 2 . Finally the 
number of obscured stars per square arcminute is 



An* = A c 



ed -n± = TT- FWHM 2 • n* 



(4) 



An* was added to the number density found on the raw 
data. For the center of Carina the correction amounts to 
2.5%; for Fornax, the most heavily crowded galaxy, this 
corrrection amounts to 13% in the very center. This ap- 
proach still neglects the further probability of having mul- 
tiple (> 2) stars overlap. 

A King ( 1966) profile was fitted to the derived density 
profile. The derived parameters are listed in Tabled It is 
important to distinguish between the King 119621 profile, 
which we also call the empirical one, and the profile as 



Table 1. Morphological and physical parameters as deter- 
mined in this work for the Carina, Sculptor and Fornax 
dSphs. 



Parameter 


Carina 


Sculptor 


Fornax 


rr \ 120001 


6h4l'34" 


lh00'28" 


2h40'4" 


S [120001 


— 50°57' 


— 33°42' 


— 34°31' 


PA f°l 


64 ± 2.5 


98 ± 2 


40 ± 5 


e 


0.32 ± 0.04 


0.2 ± 0.05 


0.3 ± 0.05 


emp r • 1 

r [arcmin] 


11.96 ± 1.5 


7.56 ±0.7 


13 ±0.15 


r e t mp [arcmin] 


22.54 ± 1.4 


40 ±4 


89 ±17 


c emp 


0.27 ±0.08 


0.72 ±0.2 


0.83 ±0.1 


/ a 2 ) theo 


1.2 


2.7 


3.8 


rf™ [arcmin] 


14 


10 


15 


r\ heo [arcmin] 


33.5 


44 


98 


^theo 


0.38 


0.64 


0.82 


[10 5 L Q ] 


(3.45 ±0.5) 


(5.6 ±0.5) 


(88 ±0.5) 


My [mag] 


-9.0 


-9.5 


-12.5 



published in 1966, called theoretical profile hereafter. The 
empirical profile is given by the simple equation 



J(r) 



1 



1 



i + (t-) 2 



i + (t:) 2 



(5) 



This empirical profile is usually parameterized in terms 
of the concentration parameter c = log(r t /r c ) as well as 
the core radius r c . On the other hand, following Binney 
and Tremaine fTSEZ}, the theoretical profile is given by 
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Fig. 2. Contour plot of the Sculptor dwarf spheroidal. The 
density levels correspond background value (dotted line), 
la above that (thin solid line), 2a, 5a, 10a and so on (thick 
solid lines) . Note the increase of ellipticity with radius and 
the potential tidal tails. 
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Fig. 3. Contour plot of the Fornax dwarf spheroidal. The 
density levels correspond to background value (dotted 
line), lcr above that (thin solid line), 2a, 5a, 10a and so 
on (thick solid lines). The inner isopleths show a steeper 
decline to the south-east than to the north-west. The lo- 
cations of the five globular clusters are indicated by their 
numbers. 

a differential equation, Poisson's equation, relating the 
gravitational potential *f?(r) and the dimensionless radius 
r* = r/ro, where ro is some typical radius called the core 
radius: 



d 
dr* 



*\2 



(r*) 



dr* 



47rGpi(r*) 2 



(6) 



The dimensionless density p* = p/ po is then given by the 
relation 



P*W = Pi 



erf 



4* 



i 



2* 
3^2 



(7) 



The profile is then parameterized by ip(0)/a 2 , po, and ro. 
Both profiles, fitted to an observed density profile, are 
indistinguishable down to a density contrast of two orders 
of magnitude from the central density. Yet the derived 
tidal radii and concentration parameters differ by up to 
a factor of f.5. We decided to use the theoretical King 
( 1966) formula, as it has a physical interpretation. This is 
opposite to the practice in IH and Majewski ( 2001, private 
communications). For comparison purposes we also quote 
the tidal radii as determined by a fit with an empirical 
King l|1962(l model, although the fit is somewhat worse 
for all three galaxies. 

For Carina the King profile fits extremely well down 
to the limits of the data (fig. 0}. The limiting radius 



r*' leo =31.8' is similar to previous studies despite the use 
of the somewhat different theoretical profile. The tidal ra- 
dius determined from the empirical profile r t = 23.6' is 
25% smaller than the IH value (rf H = 28.8). After back- 
ground subtraction negative values for the density occur 
where the profile drops to the background. As these nega- 
tive values are not depicted in a logarithmic plot, it looks 
as if the profile would reach a constant positive value, yet 
this is not the case. The seemingly constant value of the 
profile is a good measure for the standard deviation of 
the background. To illustrate this better we also show the 
absolute value of the negative values of the background 
as open circles. Although we show for the first time the 
density profile of the main sequence stars down to a den- 
sity of 4 x 10~ 3 • po (where po is the central density), we 
cannot confirm the break to a shallower slope found by IH 
and Majewski et al. (2000) at a radius of 20 arcminutes. 
This break should occur at l/25th of the central density. 
At this significance we can also rule out any unsymmetric, 
tidal-tail like form for this second component. 

For the Sculptor density profile (fig. [SJ a departure 
from the fitted profile can be seen clearly at 30' from the 
center. This break in the profile is further evidence that 
an extended stellar component exists. The extended com- 
ponent is undistinguishable from the background outside 
of about 45'; its interpretation is deferred to section 
The parameters determined from the King-fit again yield 
a smaller tidal radius and a smaller concentration param- 
eter, compared to the IH values of 76.5' and 1.12 respec- 
tively. In this case also the values from the theoretical 
profile of r t = 44' and c = 0.64 are smaller than those 
previous determinations. We do not have an explanation 
for this somewhat disturbing difference. The only sizeable 
reduction difference lies in the crowding correction, but 
we would not expect this to make such a difference. We 
choose e=0.2 in this work as this seems to be the value 
in the undisturbed center of the galaxy, whereas the IH 
value of 0.32 might be influenced by the unresolved tidal 
tails. For comparison with earlier work, the corresponding 
values derived with an ellipticity of 0.3 are: r t = 47' and 
r c = 7.6'. 

The Fornax density profile (fig. ||JJ) follows the King 
profile very well on to where it drops to the background 
noise of about 0.1 counts per arcmin 2 . The King limit- 
ing radius v t — 98' is bigger than previously determined 
(IH: r f =71', c=0.72) and the profile extends almost to the 
tidal radius derived from the King profile fit. On the other 
hand, our field of view is only marginally bigger than the 
galaxy itself. The background area is therefore not as big 
as in the two other cases, so a really very extended second 
component is not excluded by this data. 

An equally acceptable fit is also obtained with Sersic 
profiles. They are given by 



J(r) = I e e 



(8) 



We obtain an R e of 9', 8' and 16' while n is 0.6, 0.8 and 0.7 
for Carina, Sculptor and Fornax respectively. This result 
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Fig. 4. Background subtracted density profile for the 
Carina dSph. Clearly a King profile with rt = 31.8' and 
c = 0.35 fits extremely well out to the limits of the data. 
After background subtraction negative values of the den- 
sity occur which do not show up in the logarithmic plot. To 
illustrate these, we plotted their absolute values as open 
circles. 



Fig. 5. Background subtracted density profile for the 
Sculptor dSph. A King profile with rt = 70' and c = 0.76 
fits well out to 25'. Then the profile flattens and hits the 
background at about 45'. After background subtraction 
negative values of the density occur which do not show 
up in the logarithmic plot. To illustrate these, we plotted 
their absolute values as open circles. 



matches with the picture described in Caon et al. 1993 
and Jerjen et al. 2000. These authors have fitted Sersic 
profiles to big samples of elliptical galaxies and find that 
the exponent in the Sersic profile follows a relation with 
the scale length or alternatively the absolute Magnitude. 
Smaller and fainter galaxies are better fitted by a larger 
exponent, where the largest exponents they find for small 
galaxies are just smaller than 2. There is no reason to rule 
out the Sersic profile. Nevertheless the King profile will be 
used as the benchmark, because it has a physical interpre- 
tation and the required cut-off. There are numerous other 
possible profiles, but it is beyond the scope of this paper 
to explore all of them. 

To determine the total luminosities we first derived 
the luminosity function of each dwarf galaxy. To that end 
we normalized the distribution in apparent magnitude of 
all objects in a background area to the area of the whole 
catalogue (after applying the star/galaxy classifier only 
for Carina). We then subtracted this distribution func- 
tion from the distribution function of the whole object 
catalogue. All stars belonging to the dwarf being at the 
same distance this yields the luminosity function of the 
dwarf galaxy. The integral over the luminosity function 
then gives the total luminosity L^*. 

For Carina, using a distance modulus of m-M — 20.05 
(MighellEMIl), we find L*/ = (2 ± 0.5) x 10 5 L Q ). We 
excluded stars brighter than 18th magnitude to avoid con- 
tamination by bright foreground stars. Using our data 



alone we also have to impose a cut in apparent magnitude 
at 23.5 mag. To improve our estimate we extrapolated the 
luminosity functions of the three dwarfs using published 
HST data. The total luminosities we quote are therefore 
free of any assumptions regarding the shape of the lumi- 
nosity function. For Carina we used two datasets, namely 
a ground-based but deeper one from Hurley-Keller et al. 
l|TM5|l and HST data from Mighell l(TM?|) . We first nor- 
malized the Hurley-Keller luminosity function to ours and 
then extrapolated further (until approximately 27th mag- 
nitude) by normalizing the Mighell luminosity function to 
the normalized Hurley-Keller function. This is necessary, 
because the number of stars brighter than 23rd magni- 
tude is rather low in the HST data, so the intermediate 
step avoids errors due to small number statistics. The to- 
tal luminosity of Carina is then Ly* = (3.4 ± 0.5) x 10 5 
Lq. The quoted error tries to account for the unknown 
contribution from background galaxies in the HST data 
and for the very faint stars in Carina still missing in the 
integration 

For the Sculptor dwarf the same approach as in 
Carina was necessary, because the relevant HST data set 
(Monkiewicz et al. 1999 ) was originally designed to avoid 
crowding problems, so that there are again almost no stars 
brighter than 23rd magnitude. Yet the same authors find 
that the stellar content of Sculptor is similar to Fornax, 
so that we decided to interpolate the luminosity function 
with Fornax data from Buonanno et al. 119991 With a dis- 
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Fig. 6. Background subtracted density profile for the 
Fornax dSph. Clearly a King profile with r< = 98' and 
c = 0.82 fits extremely well out to the limits of the data. 
After background subtraction negative values of the den- 
sity occur which do not show up in the logarithmic plot. To 
illustrate these, we plotted their absolute values as open 
circles. The background value relies only on the two last 
points for this profile. 

tance modulus of 19.3 (IH) and a limiting magnitude of 
again about 27mag, we derive Ly 4 = (5.6 ± 0.5) x 10 5 Lq 
(from our data alone: (4.6±0.5) x 10 5 Lq). 

For Fornax the star catalogue found in Buonanno et al. 
( 1999) is sufficiently populated to be able to skip the inter- 
mediate step. Taking all stars brighter than 26th magni- 
tude this yields a total luminosity of L^' = (8.7±0.5) x 10 6 
L Q (from our data alone: (6.7±0.5) x 10 5 L Q ), where 
m - M = 20.68. 

We did not correct for stars fainter than 27th magni- 
tude; our values might therefore be considered as lower 
limits. This is consistent with all three luminosities being 
smaller than those quoted in Grebel (2000) by about 30%. 



4. Summary of observations and dicussion 

In this paper we present a new study on the morphologi- 
cal parameters of three Galactic dwarf spheroidal galaxies, 
namely Carina, Sculptor and Fornax. We show the stellar 
density distributions below the main sequence turnoffs in 
fields of 4, 16.25 and 8.5 square degrees, respectively. In 
the Sculptor dwarf we discover indications of tidal tails, 
but the quality of the current data does not allow to 
draw firm conclusions. Following the seminal paper by 
Irwin & Hatzidimitriou (1995) we expected to detect ex- 
tratidal populations of some sort in all of our three ob- 
jects. Yet we can rule out those populations for two of 



them. Together with similar findings by Odenkirchen et al. 
I|2001|) for Draco this is conclusive proof that not all dSphs 
really have extended populations on the level claimed by 
IH. Nonetheless, considering Carina alone, it remains dis- 
turbing that although our density contrast should be am- 
ply sufficient we cannot confirm the extratidal component 
found by Majewski et al. |2000| ) in their study of the RGB 
star population (see sections ri] & |3J) . This discrepancy is 
at the moment difficult to explain. It is possible that the 
extended component is so large that we do not reach the 
real background value, thereby overestimating it. Another 
possible explanation is related to the different distribution 
of differently aged populations in Carina found by vari- 
ous authors (Harbeck 2001, Monclli 2003) From our data 
alone we conclude that no direct sign of interaction with 
the tidal field of the Milky Way can be seen. This supports 
the idea that Carina might be dominated by dark matter. 

In the context of cold dark matter theory (CDM), the 
important piece of information to derive from the data 
is the mass of any dark matter halo in which the dSph 
galaxies may be imbedded. The usual approach is to de- 
termine a mass to light ratio (M/L) from kinematic and 
photometric data. There have been some doubts whether 
the kinematic approach actually determines a virialized 
mass or should rather be considered as an upper limit, as 
tidal forces might increase or even dominate the measured 
velocity dispersion (for example Klessen & Kroupa fl998fl . 
We will therefore now determine a lower limit for the M/L 
from the morphology of the dwarfs only, using a similar 
approach as Faber & Lin <|1983[> and Pryor l|1996jl . The 
King profile is strictly speaking only valid for a relaxed 
isothermal sphere with a cutoff, where no orbital motion 
of any kind is being taken into account. Nevertheless it 
can be generalized for an elliptical orbit. The perigalactic 
tidal radius can then be taken as a good approximation to 
the effective tidal radius fKing Hl?r)2*| . Using this assump- 
tion Oh, Lin & Aarseth (1992) give the following relation 
between the King tidal radius of a satellite and its mass: 



where 

m = 



[(l + e) 2 /2e] ln[(l + e)/(l-e)] + l 



(9) 



(10) 



and where e = (r a 



i)/(r a 



r P eri) and a are the 



eccentricity and the semimajor axis of the satellites orbit, 
respectively. MdSph and Mq are the total mass of the dSph 
and the mass of the Galaxy inside a sphere with radius a 
(see also Burkert 1992J). For a logarithmic potential the 
mass of the Galaxy inside a radius r is 



Mm. 



MoW = f ,1.1X10'° 



Combining eg uat ions 191 and 1 1 1 1 we have: 



M d sph = r\ x 



1.1 x 10 10 M„ 



lkpc 



x/(e) 



(11) 



(12) 
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Fig. 7. The dependence of the lower limit for M/L on 
eccentricity according to equation 1121 for orbital major 
axes of 1 and 2 times the current distance. In the up- 
per panel the normalized perigalactic distance is plotted 
against M/L. The corresponding eccentricity can be read 
off the lower panel. The chosen apogalactic distances are 
the following. Sculptor (r t = 562pc): r apo — 88kpc (solid 
line) and r apo = 176kpc (dotted). Carina(r t = 546pc): 
Ta V o — 94kpc (long dash) and r apo = 188kpc (dot - long 
dash). Fornax(r t = 3600pc): r apo — 138kpc (short dash) 
and r apo = 276kpc (dot - short dash). 



In order to derive lower limits for M/L we choose the 
lower limits for the tidal radii to be the largest radii where 
the galaxies are still obviously regular in shape. For Carina 
and Fornax this is the point where the profile drops to the 
background, so we have r™ r = 20' and r/°' = 90'. For 
Sculptor we choose the break radius in the profile as the 
lower limit, which is around rf cl = 25'. The lowest mass 
configuration is when the dSph is on a circular orbit. We 
thus derive the lowest possible limits for MdSph given in 
Table El 

Figure illustrates how MdSph depends on the peri- 
galactic distance or equivalently the eccentricity e for or- 
bital major axes a of 1.5 and 2.5 times the current distance 
d (taken from Grebel EUUUjl . 

Alternatively Table also gives a "best guess" for 
MdSph where we choose almost randomly an eccentricity 
of 0.6 and an apogalacticon of 2 times the current distance. 
For Carina and Fornax we take the tidal radius from the 
King profile fit, while for Sculptor the best guess remains 
the radius where the first asymmetries appear. 

According to Saviane et al. (2000) the main stellar 
population of Fornax is 5.4xl0 9 years old which corre- 
sponds to an M/L in the V-band of about 4.8 according 



Table 2. Mass estimates for the Carina, Sculptor and 
Fornax dSphs. The quantities labelled with mm are lower 
limits for the mass and M/L obtained by assuming a circu- 
lar orbit. The quantities labelled best give a "best guess" 
based on almost random, but conservative assumptions 
for the orbits. Also given are the heliocentric distance, the 
kinematic M/L and the heliocentric velocity. 



Parameter 


Carina 


Sculptor 


Fornax 


Distance [kpc] 


94 


88 


138 


M/L dyn 


30 


11 


4.8 


v [km/s] 


225 


100 


50 


rr n [pc] 


546 


639 


3600 


M mm [M ] 


2.0xl0 5 


3.7xl0 5 


2.7x10 


M/L" 1 ™ 


0.6 


0.7 


3 


4 est [pc] 


914 


639 


3928 


a best [kpc] 


188 


176 


276 


^best 


0.6 


0.6 


0.6 


M best [M ] 


6xl0 6 


2xl0 6 


2 x10 s 


M/L best 


17 


3.6 


23 



to Bruzual & Chariot (2000). We therefore do not need 
any dark matter halo for Fornax. By requiring M/L to 
be equal to the kinematically derived M/L we can argue 
that Fornax must almost be in its perigalacticon. Further 
approach to the Galactic center would make dark matter 
necessary to prevent the formation of tidal tails above our 
detection limit. 

The case is not nearly as clear-cut for Carina, because 
the lowest limit to M/L is so low. Yet by setting the eccen- 
tricity to 0.46 and assuming Carina to be in its apogalac- 
ticon we can easily bring M/L to be 30 as found in kine- 
matic studies. The perigalactic distance would then be 
around r per i = 35kpc. For comparison: the nearest dwarf, 
Sagittarius, is currently at a distance of 16 kpc from the 
Galactic center. Generally speaking such an eccentric orbit 
is in good agreement with the high Galactocentric velocity 
of Carina, although this radial velocity is at the same time 
obvious proof that Carina is not exactly in its apogalac- 
ticon now. Still, as stated by Burkert l|1996(l . a detection 
of tidal tails at a lower level but well inside the tidal ra- 
dius adopted here, would make it difficult to accomodate 
a dark matter halo for Carina without invoking extreme 
orbital parameters. 

In the case of Sculptor, to bring M/L up to the kine- 
matic value, the eccentricity would have to be bigger than 
0.5, which leads to r per i — 28kpc. Assuming both the kine- 
matic and the King approach to be valid we thus have to 
assume a highly eccentric orbit for Sculptor. 

The direction of the tails compared to the direction 
of the Galactic center predicts a crude direction for the 
tangential orbital motion. The tail which is on the inner 
side of the main body will precede it in direction of the 
velocity vector, while the outer tail stays behind the main 
body of the dwarf. The direction of motion of Sculptor 
should therefore be in north-south direction. As Sculptor 
is very near the south Galactic pole the sign of the proper 
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motion vector can not be inferred because it is not clear 
which tail is the inner one. The only proper motion mea- 
surement for Sculptor by Schweitzer et al. (1995) gives a 
tangential velocity vector pointing to the north-east. 

One of the propositions to solve the overabundance 
of dark satellites is to inhibit star formation in low mass 
dark matter satellites by photoionisation heating and feed- 
back. The visible dwarfs are then identified with the most 
massive (> 1O 9 M0) satellite substructures found in CDM 
simulations. The consequence is (citing Stoehr et al. 12002(1 : 
"Although there is no problem accomodating a single dis- 
rupting object like Sagittarius, it would become uncom- 
fortable if tidal stripping were detected unambiguously in 
other systems." (see also Hayashi et al. I2002|l . Counting 
Sagittarius, Ursa Minor (Martinez-Delgado et al. 121) (Jlf) 
and Sculptor, there are now three satellites of the Milky 
Way out of a total of 13 that are most likely currently 
strongly influenced by tidal forces. Note also that even 
our most massive galaxy, Fornax, has only 10 8 M©. We 
therefore conclude that our measurements do neither sup- 
port nor contradict dark matter in dSphs, but the over- 
abundance problem for Galactic satellites clearly remains 
unsolved. 
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